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THE  COMMINUTION  BEHAVIOR  OF  LIQUIDS 


1.  INTRODUCTION 

The  aim  of  this  investigation  is  to  assess  the  role  of  nonlinear  viscoelasticity  in  the 
process  of  liquid  comminution  under  explosive  stress.  Ultimately  the  goal  of  the  program  is  to 
evaluate  the  hypothesis  that  the  crucial  factor  in  determining  the  comminution  and  dispersal 
characteristics  of  liquids  under  these  conditions  is  the  nonlinear  dynamical  response  behavior  of 
the  liquids.  In  attacking  the  problem,  it  is  essential  first  to  develop  a  model  of  comminution 
that  incorporates  the  mechanism  of  liquid  response  and  failure  when  subjected  to  explosively 
generated  stresses.  The  central  theme  of  this  report  is  the  proposal  of  such  a  model,  the  ex¬ 
ploration  of  its  consequences,  and,  where  possible,  a  comparison  of  predictions  based  upon  it 
with  experimental  data.  This  model  calculation  forms  the  first  part  of  the  report.  Following 
this  is  a  brief  section  which  isolates  those  elements  of  the  proposed  mechanism  that  are  con¬ 
trolled  by  the  nonlinear  viscoelastic  characteristics  of  the  liquid.  Lastly  we  address  the  issue  of 
utilizing  computer-simulation  molecular-dynamics  (MD)  experiments  as  a  tool  for  probing  these 
nonlinear  characteristics;  we  conclude  that  MD  experiments  offer  significant  potential  for 
contributing  to  the  understanding  and  control  of  the  processes  of  explosive  liquid  comminution 
and  dispersal. 

2.  BACKGROUND 

The  experimental  study  that  provides  the  principal  impetus  and  guidance  for  addressing 
comminution  behavior  is  that  of  Gerber  and  Stuempfle.*  In  this  work  a  high-explosive  projector 
was  employed  to  produce  dissemination  of  a  liquid  cloud  around  a  preferential  vertical  axis.  A 
three-dimensional  sampling  array  was  used  to  measure  the  mass  distribution  of  the  cloud  follow¬ 
ing  its  emergence  from  the  projector.  The  mass  distribution  along  the  vertical  z-axis  was 
described  by  a  Weibull  distribution, 

f(z)  =  (/3/z)(z/z0)^  exp[-(z/zQ)^J 

such  that  f(z)dz  gives  the  mass  fraction  lost  in  the  height  range  z  to  z  +  dz  above  the 
projector.  In  the  Weibull  forir,  Zq  is  a  characteristic  length  parameter  and  /3  is  a  shape  parame¬ 
ter.  Gerber  and  Stuempfle  present  their  data  in  terms  of  four  parameters  z^q  —  the  median  of 
the  distribution  —  and  the  first  three  moments  of  the  distribution.  In  terms  of  the  parameters 
of  the  distribution 


z50  =Z0(A  2)1/0, 

and  the  mean  or  first  moment  is 

<Z>  =  zor  (l/B  )  /B  . 

Values  of  P  fall  in  the  range  from  about  0.6  to  3  with  most  of  the  more  than  30  liquids  studied 
being  described  by  p's  near  2;  <^z^>  values  varied  from  about  65cm  to  about  335cm.  Generally 
one  finds  that  larger  values  of  <z)>  are  accompanied  by  larger  values  of  p.  A  major  theme  in 
what  follows  is  the  attempt  to  understand  those  liquid  properties  that  are  operative  in  the 


*  Gerber,  B.V.  and  Stuempfle,  A.D.  TSD  004.7-1 1 /Vol.  1.  A  New  Experimental  Technique  for 
Studying  the  Explosive  Comminution  of  Liquids.  1976  Army  Science  Conference 
Proceedings.  Vol.  I.  September  1976.  Unclassified  Report. 
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determination  of  the  average  height  <z>.  Exploring  the  details  of  the  mass  distribution  pattern 
is  beyond  the  scope  of  this  work. 


A.  PRELIMINARY  CONSIDERATIONS 

In  this  section  we  give  a  crude  description  of  the  liquid  comminution  process.  The 
principal  result  of  the  section  is  an  equation  giving  the  average  height  <z>  in  terms  of  the  size 
of  the  particles  or  droplets  that  compose  the  cloud. 

We  consider  the  effect  of  the  explosion  on  the  liquid  to  be  roughly  separable  into  two 
major  components: 

(a)  The  explosive  release  of  energy  generates  a  shock  wave  that  in  passing  through 
the  liquid  produces  large  tensile,  compressional,  and  shear  stresses.  These  cause  local  internal 
fissures  within  the  liquid  so  that  it  comminutes  into  many  small  liquid  fragments. 

(b)  The  explosion  produces  a  rapidly  expanding  gas  that  intrudes  itself  into  the 
stress-produced  fissures,  thereby  preventing  their  healing,  and  causing  the  liquid  fragments  to 
be  expelled  from  the  projector  with  a  vertical  velocity  Vq. 

The  dependence  of  the  initial  velocitv  on  liquid  parameters  can  be  estimated  by  the  fol¬ 
lowing  sort  of  rough  argument.  We  shall  assume  that  relatively  little  of  the  energy  of  explosion 
is  used  in  the  formation  of  the  droplets  so  that  most  of  the  energy  E  is  transferred  to  the 
kinetic  energy  of  the  droplets.  The  initial  velocity  imparted  to  the  droplets  then  becomes 


vQ  =  (2e/m)1/2*  p"1/2  (1) 

where  m  is  the  mass  of  a  droplet,  G  is  the  kinetic  energv  per  droplet,  and  p  is  the  liquid 
densitv.  Assuming  that  E  is  the  same  for  each  explosion,  it  follows  that  the  initial  velocity  of 
the  droplets  is  dependent  onlv  on  the  density  of  the  liquid. 

The  dynamical  problem  for  the  projected  droplets  can  now  be  worked  out 
straightforwardly,  at  least  in  lowest  order  approximation.  We  assume  that  the  droplets  are 
spherical  and  that  the  size  distribution  is  characterized  by  a  function  f(a)  which  gives  the  proba¬ 
bility  of  a  droplet's  radius  falling  in  the  range  a  to  a+da.  We  take  the  forces  acting  on  a  droplet 
of  radius  a  to  be 


and 


(a)  gravity:  -i-Tra^pg, 

(b)  air  resistance:  -6xyav. 


In  these  expressions  g  is  the  acceleration  due  to  gravity,  y  is  the  viscosity  of  the  medium  (airl 
into  which  the  droplet  cloud  is  projected,  and  v  is  the  instantaneous  droplet  velocitv  in  the 
vertical  direction.  Newton's  law  then  takes  the  form 


dv/dt  +  v/T  =  -g, 


(2) 


where  for  a  droplet  of  radius  a 


T  =  2pa2/9ry 


ft 


(31 
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of  radius  a  to  be 


and 


(a)  gravity:  -i-Tra^pg, 

(b)  air  resistance:  -6xyav. 


In  these  expressions  g  is  the  acceleration  due  to  gravity,  y  is  the  viscosity  of  the  medium  (airl 
into  which  the  droplet  cloud  is  projected,  and  v  is  the  instantaneous  droplet  velocitv  in  the 
vertical  direction.  Newton's  law  then  takes  the  form 


dv/dt  +  v/T  =  -g, 


(2) 


where  for  a  droplet  of  radius  a 


T  =  2pa2/9ry 


ft 


(31 


The  solution  of  (2)  is  just 


(4) 


V  =  v0  exp(-t/T)-gT  [1  -  exp  (-t/x)]  . 

Integrating  this  again  gives  the  equation  of  motion 

z(t)  =  v0x [1-exp  (-t/x)  +  gx2  1  -  -£■  -exp  (-t/x)].  (5) 

The  height  to  which  a  droplet  rises,  z,  can  then  be  found  by  setting  v=0  in  (4)  and  solving  for  the 
time: 


t0  =  x  i,n(l  +  vQ/gX)  (6) 

Putting  z  =  z  (t  )  gives 

o 

2 

z=v0x-gx  UnO  +  v0/gx).  (7) 

<1 

For  typical  values  of  »?,  p,  and  a  (4=  0.18  mP,  Psil.3  g/cm  ,  acs:25  pm),T—  .01  sec.  Using  vQss  3 
x  104  cm/sec  and  g  =  980  cm/sec''  one  gets  t  —  .08  sec,  which  is  of  the  correct  order  of 
magnitude  to  what  is  experimentally  observed.  Moreover  vQ»gr  so  that  approximately 

z  -  vqX.  (8) 

Equation  (8)  is  accurate  to  better  than  a  percent  or  so  for  the  full  range  of  liquids  and  experi¬ 
mental  conditions  studied  by  Gerber  and  Stuempfle.  Using  equation  (3)  and  averaging  gives 

<z>  =  (2v0/9n)p<a2>,  (9) 

where  ‘(a2)  is  the  mean  square  droplet  radius. 


4.  A  PHYSICAL  MODEL  FOR  COMMINUTION 

The  central  issue  is  clear  from  equation  (9):  What  are  the  processes  -  and  thus  the 
liquid  characteristics  -  that  determine  the  mean  square  droplet  size?  We  assume  that  the 
mechanism  of  droplet  formation  is  mainly  a  mechanical  (as  opposed  to  a  thermal)  one.  At  any 
instant  even  an  equilibrium  liquid  exhibits  local  deviations  from  homogeneity  and  isotropy.  One 
characterization  of  the  deviations  from  homogeneity  is  the  distribution  of  fluctuations  in  the 
density.  From  classical  statistical  mechanics  the  probability  of  occurrence  of  a  density 
fluctuation  in  the  range  AptoAp+  dp  is 

P(Ap)dP  =  (2tt<  Ap2^  exp[-ApZ/Z<  Ap2>]  dp,  (10) 

,  2 

where  <  Ap  >  is  the  mean  square  density  fluctuation  given  by 

<  Apz>  =  p2  (kT/V)/K,  (11) 


where  T  and  p  are  the  temperature  and  average  density  of  the  system,  V  is  the  volume  of  the 
region  under  consideration,  k  is  Boltzmann's  constant,  and  K  is  the  equilibrium  bulk  modulus  of 
the  material. 

We  suggest  that  under  the  influence  of  the  large  explosively  generated  stresses,  these 
local  inhomogeneities  will  behave  as  structural  flaws  in  the  material:  shear  stresses  tend  to 
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elongate  the  inhomogeneities  by  "organizing"  the  liquid  structure*  while  the  tensile  stresses 
cause  the  flaws  to  propagate  at  the  stress  intensification  points.  The  situation  is  depicted 
pictorially  in  figure  1.  However,  in  a  liquid,  even  at  the  elevated  pressures  produced  by  the 
explosive  shock,  one  expects  that  the  rate  at  which  the  flaws  heal  and  the  rate  at  which  thev 
propagate  via  growth  at  their  tips  should  be  comparable;  consequently,  one  is  not  dealing  with  a 
classical  brittle  fracture  situation,  but  rather  one  involving  relatively  little  growth.  Thus  the 
maximum  flaw  dimension  should  be  comparable  with  the  size  of  the  particles  that  result  from 
the  explosive  comminution. 

It  is  the  processes  described  in  this  last  paragraph  that  will  depend  sensitively  on  the 
nonlinear  viscoelastic  response  characteristics  of  the  liquid  system.  Because  the  large  stresses 
to  which  the  material  is  subjected  are  "switched  on"  essentially  instantaneously,  the  nature  of 
the  behavior  will  be  governed  in  large  measure  by  the  characteristic  time  scales  of  the  shear 
and  structural  responses  of  the  liquid,  i.e.,  the  shear  and  structural  relaxation  times.  The 
specific  picture  that  we  have  described  above  is  an  oversimplified  one  in  which  we  have  as¬ 
sumed  that  the  .hear  relaxation  time  is  negligibly  small  in  comparison  with  the  structural 
relaxation  time.  Only  in  this  limiting  situation  is  one  able  to  portray  the  flaw  development  and 
liquid  fracture  processes  without  reference  to  the  specific  viscoelastic  and  strength  parameters 
that  describe  the  material.  Generally  speaking  these  extreme  conditions  do  not  obtain;  the 
shear  relaxation  time,  while  smaller  than  the  structural  time,  is  of  the  same  order  of  magni¬ 
tude.  Consequently,  the  sequential  nature  of  the  mechanisms  that  operate  in  comminution  that 
was  assumed  above  —  the  development  of  large  shear  strains,  followed  by  structural  reorganiza¬ 
tion  (flaw  development),  followed  by  liquid  fracture  from  tensile  stresses  —  is  not  nearly  so 
precisely  defineable.  The  various  mechanisms  do  not  proceed  strictly  serially,  but  rather  couple 
in  a  complicated  interactive  and,  perhaps,  competitive  fashion.  Depending  upon  the  relative 
time  scales  and  specific  forms  of  the  various  response  mechanisms,  this  coupling  can  either 
reinforce  the  liquid  comminution  process  (leading  to  smaller  particle  sizes)  or  interfere  with  it 
^resulting  in  larger  particle  sizes).  Nevertheless,  the  elementary  picture  that  we  are  examining 
is  a  useful  one  for  it  provides  a  general  framework  in  which  to  analyze  and  compare  the 
behavior  of  different  liquids  acted  on  by  explosive  stresses,  as  well  as  highlighting  the  crucial 
liquid  properties  that  govern  (at  least  in  lowest  order  approximation)  the  cloud  dispersal  charac¬ 
teristics. 

In  order  to  estimate  the  droplet  size,  we  can  employ  the  following  sort  of  elementary 
notions.  We  consider  that  in  order  to  generate  droplets  of  characteristic  dimension  =  a,  a  region 
of  the  liquid  of  volume  Vssa'5  must  contain  a  flaw  whose  volume  exceeds  some  critical  volume 
Vc<  The  action  of  the  explosive  shear  stresses  is  to  elongate  and  shape  the  flaw  [into  roughly  a 
"pillow"  shape  (see  figure  2)]  so  that  it  can  be  described  hv  dimensions  a,  b,  and  c  (a>b>c).  The 
minimum  flaw  volume  V  =  abc,  that  will  lead  to  the  production  of  droplets  with  dimensions  on 
the  order  of  "a"  will  result  if  one  has  a  local  density  fluctuation 

Ap  >_  pVc/v  =  A  pc  (12) 

in  the  volume  V.  The  liquid  will  comminute  into  particles  of  this  size  only  if  a  fluctuation  of 
such  magnitude  occurs  with  significant  probability.  The  condition  for  this  to  obtain  is 

Apc«<  Ap2>1/2:  (13) 
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Figure  l.A.  A  Schematic  Representation  of  Liquid  Flaws  Elongated  under  the  Action  of  Large 

Shear  Stresses 


Figure  1. 


.  Liquid  Flaws  That  Have  Grown  and  Intersected  to  Form  Comminution  Droplets 
under  the  Action  of  Tensile  Stresses. 


that  is,  the  critical  density  fluctuation  must  be  comparable  with  the  PMS  density  fluctuation. 
Substituting  equations  (11)  and  (12)  in  equation  (13)  and  using  Vssa3  and  V  =  abc  leads  to  an 
approximate  expression  for  the  value  of  a: 


a  -  (be)2  K/kT. 

Recalling  equation  (9),  this  leads  to  an  equation  for  the  mean  height  of  the  droplets: 


(14) 


<z>  =  (2/9n)  vQp  (be)4  (K/kT)2.  (15) 

In  order  that  the  flaws  do  not  heal  themselves  too  rapidly,  it  is  essential  that  the  smallest  flaw 
dimension  c  be  at  least  on  the  order  of  a  molecular  separation  distance;  with  this  in  mind  we 
have 


-1/3 

c  o  P  (16) 

-1/2 

and  recalling  that  v0  »  p  '  ,  we  may  collect  together  all  the  density-dependent  quantities  on 
the  right  hand  side  of  equation  (15).  The  resulting  expression  takes  the  form 


<z>  oc  p-5/6  (17) 


5.  COMPARISON  WITH  EXPERIMENT 

Equation  (17)  is  in  a  form  that  can  be  compared  with  experimental  results.  The 
comparison  is  shown  in  figure  3  for  six  liquids  at  20°C.  Many  of  the  data  for  this  figure  were 
taken  from  the  work  of  Gerber  and  Stuempfle*  and  are  presented  in  the  table  which  will 
follow.  As  is  apparent  there  is  fairlv  good  agreement  between  these  data  and  the  prediction  of 
equation  (17)  (the  slope  of  the  line  shown  in  figure  3  is  65.5  g^ cm-^' 2GPa~2).  This  offers 
fairly  strong  support  for  the  essential  correctness  of  the  view  of  the  comminution  mechanism 
that  was  described  above.  It  is  also  clear  that  the  modifications  of  this  overly  simplified  view 
that  are  required  if  one  correctly  treats  the  dynamical  aspects  of  the  liquid  response  are  not 
understood  in  any  kind  of  quantitative  fashion.  Some  of  the  molecular  dynamics  results 
(presented  in  the  appendix)  give  a  fair  indication  of  the  type  of  phenomena  that  must  be 
considered  —  especially  in  regard  to  liquid  structural  changes  occasioned  by  large  shear  stress  — 
but  efforts  to  apply  these  results  quantitatively  to  comminution  are  at  this  point  premature. 

Of  the  six  liquids  which  were  considered  above,  except  for  carbon  tetrachloride, 
comminution  data  have  also  been  obtained  at  130°C.  The  variation  that  one  would  expect  from 
our  model  is  given  by  equation  (15);  the  major  effect  appears  to  be  the  factor  (1/kT)  .  On  the 
basis  of  this  one  would  expect  the  ratio  of  the  mean  height  at  20°C  to  that  at  130°C  would  be 
about  (403/293)2%1.9.  The  experimentally  determined  ratios  for  the  liquids  tha*  we  have  con¬ 
sidered  fall  into  the  range  of  about  1.0  to  1.5. 
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The  Mean  Height  to  Which  the  Droplets  Rise  Multiplied  hv  Density  to  the  5/6 
Power  Plotted  Versus  the  Square  of  the  Bulk  Modulus.  The  Solid  Line  Is  Drawn 
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Table.  Data  Used  in  Constructing  Figure  3.  All  Data  Are  Given  for  a  Temperature  of  20 °C. 


Liquid 

Density 

<z> 

Modulus 

Tetrachloroethvlene  (  o  ) 

(g/cm3) 

1.62 

(cm) 

64 

(GPa)a 

1.8 

Carbon  tectrachloride  (  •  ) 

1/59 

96 

1.39 

2,2-dichloroethanol  (A) 

1.40 

232 

1.91  ; 

1 ,2,3-trichloropropane  (A) 

1.61 

182 

2.03 

Water  (□) 

1.00 

336 

2.24 

Tetrabromoethane  (•) 

2.91 

201 

2.81 

(a)  lGPa  =  109N/m2  =  1010dyne/cm2  =  145.5  kpsi 

This  cannot  be  accounted  for  by  considering  the  change  with  temperature  of  the  density  and  the 
modulus.  It  probably  reflects  the  increasingly  important  role  of  the  viscoelastic  responses,  for 
which  the  relaxation  times  will  be  significantly  shorter  at  the  higher  temperature. 


6.  THE  EFFECT  OF  POLYMERIC  ADDITIVES 

Gerber  and  Stuempfle*  also  report  measurements  at  20°C  on  two  samples  of 
tetrachloroethvlene  to  one  of  which  was  added  1  %  polvisobutylene  (PIB),  while  the  other 
contained  1%  polystyrene  (PS).  In  the  former  case  <z>was  increased  by  nearly  a  factor  of  four 
while  in  the  latter  case  the  increase  was  less  than  50%.  Phrased  another  way,  the  effect  of  PIB 
is  nearly  to  double  the  average  droplet  radius,  whereas  PS  produces  only  about  a  20%  increase  in 
droplet  size.  The  model  of  comminution  that  we  have  considered  would  lead  one  to  expect  am 
increase  in  droplet  size  on  the  order  of  10%  or  less.  Clearly  some  other  mechanism  is  operative 
in  determining  the  particle  size.  We  suggest  that  the  effect  of  the  polymeric  molecules  is  to 
act  as  agents  that  terminate  the  flaws  bv  inhibiting  their  growth  under  the  applied  stresses. 
The  PIB  molecules,  which  were  of  substantially  greater  average  chain  length,  are  much  more 
effective  in  this  regard  than  the  (relatively  shorter)  PS  molecules,  and  thus  lead  to  the 
formation  of  rather  large  droplets. 


7.  SUMMARY  AND  CONCLUSIONS 

In  this  report  we  have  presented  a  physical  model  of  the  explosive  comminution  process 
in  liquids.  The  key  feature  of  the  model  involves  the  development  of  flaws  in  the  liquid  under 
the  action  of  large  explosively  generated  stresses.  It  is  at  these  flaws  that  the  material  fails 
via  a  process  that  we  have  termed  "liquid  fracture."  Guided  bv  previous  molecular  dynamics 
(MD)  studies  that  we  have  carried  out  on  the  nonlinear  response  character  of  liquids  (a  partial 
summary  of  these  is  given  in  the  Appendix),  we  have  proposed  a  highlv  simplified  view  of  this 
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process.  Even  this  simplified  picture  is  relatively  successful  in  yielding  semiquantitative  agree¬ 
ment  with  experimental  data. 

In  order  to  develop  a  truly  quantitative  picture,  it  is  our  belief  that  the  dynamical 
aspects  of  the  liquid  failure  mechanisms  must  be  more  fully  understood.  A  thoroughgoing 
molecular  dynamics  investigation  of  the  time-dependent  structural  response  of  liquids  to  high- 
amplitude  shear  and  tensile  stresses  would  prove  most  useful  in  elucidating  the  essential 
features  of  these  nonlinear  viscoelastic  dynamics  that  operate  in  explosive  comminution.  We 
are  at  present  examining  the  possibilities  of  initiating  such  an  effort. 
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APPENDIX 


MOLECULAR  DYNAMICS 


Not  only  in  the  process  of  comminution  but  also  in  a  rather  wide  variety  of  technological 
applications,  for  example  elastohydrodynamic  lubrication,  shock  loading,  solid  propellant 
failures,  etc.,  liquids  and  amorphous  solids  are  subjected  to  extremely  large  shearing  forces. 
Because  the  microscopic  mechanisms  of  importance  are  not  well  understood,  solutions  to 
problems  involving  nonlinear  behavior  are  quite  limited. 

For  example  in  order  to  predict  shear  failure  (or  tearing)  of  a  liquid  under  the  forces  of 
an  explosion  or  shock  wave,  one  must  be  able  to  compute  the  time  evolution  of  the  stress  in  a 
fluid  element  when  it  encounters  a  simultaneously  applied  time  varying  shear  rate  and  a  time 
varying  pressure.  The  response  of  the  system  under  these  conditions  can,  for  low  shearing  rates, 
be  parameterized  in  terms  of  the  shear  viscosity,  q,  a  shear  rigidity  modulus,  Goo,  and  a 
distribution  of  relaxation  times.  At  the  rather  large  shearing  rates  generally  encountered,  this 
relatively  simple  description  must  be  modified  to  incorporate  the  nonlinear  character  of  the 
dynamical  stress  response.  Some  authors*  suppose  that  this  can  be  done  by  appending  the  idea 
of  a  limiting  shear  stress  to  the  low  amplitude,  i.e.,  linear,  shear  response.  In  order  to  test  this 
hypothesis  and,  more  generally,  to  provide  a  framework  on  which  to  base  analytical  repre¬ 
sentations  of  nonlinear  viscoelastic  behavior,  we  have  used  Molecular  Dynamics  (MD)  to  inves¬ 
tigate  a  model  liquid  system  under  conditions  approximating  those  found  in  at  high  shear  rates. 
Even  though  the  model  system  is  particularly  simple,  we  have  been  able  to  reproduce  the  major 
features  of  shear  response  of  real  liquids.  Because  an  MD  experiment  provides  one  with  a 
complete  microscopic  record  of  the  system's  evolution  in  time,  one  can  perform  "experiments" 
and  "measure"  aspects  of  its  behavior  that  are  not  accessible  in  conventional  experiments,  but 
which  provide  useful  insights  into  the  origin  of  the  properties  of  liquids  under  extreme  stress 
conditions. 

The  MD  calculation  procedure  observes  a  portion  of  liquid  as  it  is  sheared.  Although  a 
simultaneously  applied  time  dependent  pressure  and  shear  rate  would  more  closely  represent  the 
conditions  of  a  liquid  undergoing  explosive  comminution,  we  have  confined  our  attention  up  to 
now  to  considering  the  response  of  a  model  fluid  to  a  time  varying  shear  rate  at  a  series  of 
densities  (pressures). 

THE  MODEL 

MD  is  a  computer  simulation  technique  for  probing  the  microscopic  behavior  of  a  system 
that  allows  one  to  account  for  the  many-body  nature  of  the  molecular  interactions.  It  involves 
following  some  representative  number  of  molecules,  N,  by  solving  the  classical  equations  of 
motion  on  a  high  speed  computer  to  obtain  their  trajectories  in  time  steps  of  typically  0.01  ps. 
Since  only  rather  a  limited  number  of  molecules  is  treated  (generally  less  than  a  thousand),  the 
molecules  are  confined  in  a  cubic  box  which  is  surrounded  by  images  of  itself  to  avoid  severe 
boundary  effects.  In  the  studies  described  below  the  artificialities  introduced  by  this  procedure 
are  rather  insignificant.  From  the  results  of  an  MD  "experiment"  —  essentially  a  record  of  the 
system's  path  through  phase  space  —  one  can  determine  any  measurable  property  of  the  system, 
as  well  as  many  that  are  experimentally  inaccessible. 
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The  model  system  we  consider  here  is  an  assembly  of  108  Lennard-Jones  (LJ) 
particles;  these  molecules  interact  through  pairwise  additive  forces  described  by  the  LJ  6-12 
potential, 


<t>(r)  =  -4u 


(1) 


where  r  is  the  separation  of  an  interacting  pair,  and  u  and  a  are  constants  with  dimensions  of 
energy  and  distance,  respectively.  All  the  computed  quantities  are  in  so-called  LJ  reduced  units 
which  are  given  in  terms  of  u,  a,  and  m,  the  mass  of  an  individual  particle.  Although  no  real 
liquid  can  be  characterized  by  a  potential  of  the  LJ  form  (only  the  noble  gas  fluids  Ar,  Kr,  etc., 
are  quantitatively  represented  in  this  way)  their  structural  and  dynamical  behavior  can  be 
expected  to  be  (and  we  have  found  actually  are)  qualitatively  similiar  to  that  of  the  model 
fluid.  As  a  result  one  can  expect  that  the  microscopic  insights  derived  from  MD  experiments 
will  serve  as  useful  guidelines  for  understanding  and  modeling  the  behavior  of  real  liquids  and 
perhaps  for  the  "molecular  engineering"  of  improved  comminution  properties  of  liquids.  For  the 
purpose  of  comparison  with  experiment  frequent  conversion  from  reduced  to  "real"  units  will  be 
made,  using  the  parameters  for  Ar  which  are  given  in  table  A-l. 


Table  A-l.  A  Summary  of  the  Reduced  Units  in  Terms  of  Which  All  Quantities  Are  Given. 

Boltzmann's  constant  is  denoted  by  kg. 


Ouantity 

Reduced  Unit 

SI  Unit  for  Ar 

distance 

a 

0.3405  nm 

mass 

m 

6.64  x  10~2^  kg 

energy 

u 

1.65  x  10~21  J 

time 

a(m/ u)1/2 

2.16  ps 

density 

a'3 

42.1  kg  mol/m3 

temperature 

ll/kg 

119.8  K 

pressure,  stress 

u/a3 

41.8  MPa 

modulus 

u/a3 

41.8  MPa 

viscosity 

(mu)M/a2 

9.03  x  10"5  Pa  s 

thermal  conductivity 

k.g  (m/u)  ""^a-2 

1.88  x  10"2  J  m'1  K"1  s'1 

Two  methods  of  shearing  the  box  of  molecules  can  be  used.  The  results  of  each  metho/1 
agree  within  statistical  error;  however,  each  method  is  well  suited  to  examine  a  particular 
aspect  of  the  investigated  phenomenon. 


The  first  method  attempts  to  mimic  a  situation  in  which  the  fluid  forms  a  thin  film. 
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The  first  method  attempts  to  mimic  a  situation  in  which  the  fluid  forms  a  thin  film. 
This  thin  film  is  achieved  in  the  model  system  by  employing  periodicity  in  the  x  and  y  directions 
onlv;  the  external  shear  forces  are  applied  by  two  fluid  layers  translating  with  the  desired  wall 
velocities,  Uj  and  U^.  Molecules  from  the  three  regions  (i.e.,  the  two  wall  regions  and  the 
liquid  itself)  are  not  allowed  to  mix  and  are  kept  from  doing  so  by  reflection  boundary 
conditions  in  the  z  direction,  "^he  fluid  walls  (FW)  are  maintained  at  a  preset  temperature  bv 
adjusting  their  velocities  at  each  time  step  of  the  integration  scheme. 


Append  i  x 


18 


The  second  method  is  a  modification  of  the  Homogeneous  Shear  technique*  (MHS),  which 
is  better  suited  to  study  the  viscoelastic  behavior  of  the  fluid  because  it  employs  periodicity  in 
all  directions  so  that  a  bulk  material  is  sheared.  It  also  enables  a  desired  strain  rate  to  be 
instantaneously  achieved  throughout  the  liquid.  The  molecular  trajectories  are  disturbed  from 
those  governed  by  equilibrium  fluctuations  by  imposing  a  shear  strain  rate,  e,  on  the  system. 
This  is  accomplished  by  displacing  (after,  each  time  step,  At)  the  x-coordinate  of  the  i'th  mole¬ 
cule  in  the  MD  cell  by  an  amount  6x^  =  ez^At,  where  z^  is  the  z  coordinate  of  particle  i  and  At 
is  the  length  of  the  time  step,  for  as  long  as  the  shearing  is  required.  The  results  of  this 
technique,  in  particular  those  relating  to  viscoelasticity  (derived  using  a  time  dependent  strain 
rate),  complement  those  of  the  FW  method. 

We  have  used  non-equilibrium  MD  experiments  to  investigate  the  linear  and  nonlinear 
shear  response  of  the  LJ  system  under  a  variety  of  thermodynamic  conditions.  The  tempera¬ 
ture,  T,  and  number  density,  p,  initially  studied  were  close  to  the  normal  freezing  temperature, 
T  >  of  the  LJ  liquid,  that  is  0.722  k^T/u  and  0.8442  Na^/V,  respectively.  Densifications  of  10%, 
20%,  and  50%  relative  to  the  starting  density  were  achieved  by  suddenly  (within  one  time  step) 
compressing  system  at  constant  temperature.  The  densifications  were  achieved  so  rapidly 
that  the  systems  were  kinetically  prevented  from  crystallizing. 


MI.)  "M  EASUh  EVENT"  OF  SHEAR  VISCOSITY  OF  HIGH  STRAIN  RATES 


We  first  rr nsider  the  dependence  of  shear  viscosity  on  shear  rate  in  the  nonlinear  re¬ 
gion.  In  the  Homogenous  Shear  method  a  uniform  x  velocity  profile  was  set  up  instantaneously 
whereas  a  period  of  several  time  units  was  required  for  the  fluid  walls  to  drive  the  central 
region  into  this  state.  An  x  velocity  gradient  was  achieved  by  maintaining  the  upper  and  lower 
walls  at  velocities  of  U  and  -U,  respectively.  An  example  of  a  vx(z)  so  derived  is  shown  in 
A-l.  The  value  of  e  corresponds  to  an  extremely  large  shear  rate  —  on  the  order  of 
10  s  ,  depending  somewhat  upon  the  particular  values  of  u,  a,  and  m  that  are  chosen.  Such 
large  values  of  £  were  chosen  since  it  is  desirable  to  explore  the  system's  behavior  at  shear 
rates  on  the  order  of  the  reciprocal  of  the  viscoelastic  relaxation  times.  For  the  I.J  fluids  these 
are  approximately  unity  (about  10  ^  s).  The  shear  stress  needed  to  define  a  shear  viscosity  is 
readily  derived  from  the  velocities  and  positions  of  the  molecules.  "The  aB  component  of  the 
stress  tensor  is  defined  below, 
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where  V  is  the  volume  of  the  MD  cell,  v®  is  the  a  velocity  component  of  molecule  i  relative  to 
the  average  velocity  and  o  -  is  the  o  component  of  the  separation  between  molecules  i  and  j. 
The  subscripts  for  the  shear  stress  component,  <rxz,  are  dropped  in  further  discussion. 

At  low  shear  rates  plots  of  shear  stress  versus  strain  rate,  t  ,  are  linear;  however,  at 
high  stress  levels  (~Goe,  / 1 00)  the  shear  stress  increases  less  rapidly  than  (  .  In  other  words  the 
shear  viscosity  (=  oft  )  decreases  from  the  equilibrium  lor  zero  shear  rate)  value  of  tj0  with 
increasing  shear  rate.  The  MHS  method  was  used  to  obtain  y  versus  (  for  a  variety  of 
densities.  The  forms  of  v/v0  versus  stress  in  figure  A-2  resemble  closelv  those  obtained  from 
twin-disk  experiments  for  the  fluids  I>P4E,  and  the  lubricants  LA3/1271  and  Oxilube  85/140  after 
making  the  appropriate  conversions  to  real  units.  The  ratio  V/Vg  descends  markedlv  from  unity 
in  the  region  —  50MPa.  The  decrease  in  tj  with  t  is  not  due  to  thermal  heating  because  the 
calculations  were  conducted  isothermallv. 


*Hvyes.  1).  M.  cl  ah.  op.  cit. 
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I  igure  A-l.  The  z  Dependence  of  the  Average  x  Velocity  for  the  p  =  0.8442  and  i  =  0.337  3 
State  Using  the  FW  Model.  The  z  coordinate  system  is  measured  from  the 
center  of  the  film. 
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Figure  A-2. 


The  Nonlinear  Behavior  of:  A  3P  L63/1271  (P  =  0.8  GPa.  T  =  27 °C,  U  =  1.12 
m/s),  O  5P4E  (P  =  0.45  GPa,  T  =  27 °C,  U  =  0.6  m/s)  and  vOxiluhe  85/140  (P  = 
1.2  GPa,  T  =30°C,  U  =  2.2  m/s)  from  Twin-Disk  Experiments.  The  MHS 
Model:  A  p  =  0.8442,  ▼  p  =  0.92862,  •  p  =  1.01304  and  open  square  p  = 
1.2663. 
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CALCULATION  OF  SHEAR  INDUCED  STRUCTURAL  CHANGES 


The  results  in  the  previous  section  suggest  that  the  model  and  real  liquids  might  share  a 
i  umni.j.i  mechanism  of  "failure"  at  high  levels  of  shear  stress.  With  the  support  of  the  previous 
results  we  now  consider  associated  structural  changes  produced  by  the  large  shear  rates  in  the 
fluids  in  order  to  better  understand  the  failure  mechanism. 

A  directional  number  density  function,  n(a),  has  been  calculated  for  the  purpose  of 
investigating  structural  anisotropy.  The  function  n(x),  for  example,  is  the  average  number  of 
mob  rules  from  a  reference  particle  in  the  distance  range  |x|  to  |x|  +  Ax,  where  Ax  was  taken  as 
(  COO,  selected  from  the  N-l  particles  in  the  re-orientated  MD  cell.  Figure  A-3  shows  the  n(tr) 
f  the  unsheared  liquid  subtracted  from  that  of  the  sheared  system,  5n(a).  Although  the  maxi¬ 
mum  displacement  is  10To  at  most,  there  is  evidence  of  a  tendency  of  the  molecules  to  re-order 
r h. uiisel ves  into  xv  layers  (layers  normal  to  the  velocity  gradient).  The  xz  shearing  increases 
tm-  probability  of  the  molecules  in  an  xz  plane  being  found  at  multiples  of  the  intermolecular 
iiumeter  from  ra»  b  other  in  the  z  direction.  It  is  energetically  favorable  for  these  layers  to  be 
staggered  in  the  y  direction.  Although  not  so  well  defined,  these  trends  are  present  at  the 
lower  densities.  The  (grossly  simplified)  picture  that  emerges  is  one  of  corrugated  sheets  of 
atom.',  sliding  past  one  another.  A  pictorial  representation  of  this  suggested  structure  is  given  in 
figure  A  1  Th  >  structural  reorganization  results  from  the  action  of  the  imposed  velocity 
profile  which,  at  the  shear  rates  considered,  dominates  the  molecular  motion.  Major  move¬ 
ments  are  confined  to  the  line  of  shear. 

Another  aspect  of  the  shear  induced  structural  distortion  is  presented  by  the  time 
dependent  pair  radial  distribution  function,  g  (r,t),  which  is  conveniently  decomposed  into  the 
following  component  configurational  averages: 

g(£,t)  =  g0(r)  +  ( x  2 / r 2 1 uxx ( r , t )  +  ( v2 / r 2 ) u v y ( r  ,  t  ) 

+  ( zZ / r2 ) uzz ( r  ,  t )  +  ( xz / r2 ) uxz ( r , t ) 


=  .  .  .  ,  (  3  ) 

where  gp(r)  is  the  pair  radial  distribution  function  of  the  unsheared  medium. 

2 

If  f (r,t)  is  the  average  of  (  °ij^ij^rij  ^  t^le  raf^a'  element  r— *r  +  dr  about  a 
molecule  i  at  time  t,  then 

g«a(r,t)  =  15  foa(r,t)  /  (4  v  P  r2f1r)> 

and 

ga)3(r,t)  =  15  f^^(r,t)  /  (4?rpr2dr),  (4) 

-  here 

The  fafl(r,t)  have  been  calculated  at  various  times  during  the  segment.  The  difference 
between  the  g  (r,tl  of  the  steady  state  sheared  system  and  the  unsheared  svstem  at  the  same 
density,  6gQr,fr,t),  is  shown  in  figure  A-5.  It  reveals  that  shearing  produces  a  net  movement  of 
particles  in  each  coordination  shell  towards  the  origin  molecule. 

The  -.hear  rate  dependent  gxz(r,t)  at  steady  state  of  the  p  =  0.84  svstem  are  compared 
with  the  spherically  averaged  p(r,t)  from  these  samples  in  figure  A-6.  Thev  give  convincing 
vidence  of  angular  distortion  in  each  coordination  shell  such  that  in  the  positive  xz  quadrants 
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Figure  A-4.  A  Pictorial  Representation  of  the  Structural  Changes  That  Take  Place  on  Going 
from  an  Unsheared  (a)  to  a  Sheared  (b)  State.  The  arrows  denote  the  line  and 
plane  of  shear. 
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the  inside  of  ea<  li  shell  is  on  average  depleted  of  particles  when  compared  with  the  unshoared 
liquid.  In  contrast,  the  outside  of  each  shell  has  an  excess  of  molecules  in  the  positive  xz  quad- 
mis,  when  compared  with  the  unsheared  medium.  The  opposite  changes  take  place  in  the 

oeoative  i  ilia  dr  a  n  ts. 


n-  I'  !■'(  "i  OF  STRUCTURAL  REORGANIZATION  ON  MOLECULAR  MOBILITY 

It  is  to  tie  expected  that  the  effect  of  the  shear  should  not  be  confined  to  structural 
•  t.  iiges  hut  must  also  alt<  r  the  self-dif fusional  dynamics  of  the  molecules.  To  assess  this,  the 
x,  v,  and  z  components  of  the  mean  square  displacements,  <(r  (t))>,  which  exclude  the  imposed 
flow,  have  been  calculated,  i.e., 


y,y t>>  =  X>  if  t  v  ■<  «•)  dt1  -  «  (o)  ]  2>  /  n. 

o  1  t 
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tionu!  diffusion  coefficients, 
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.  s>  lection  of  the  states  are  presented  table  A-2.  The  value  of  Da  ■-  0.032  for  the  unsheared 
li  ji.ii.  near  the  triple  point  compares  favorably  with  the  equilibrium  value  of  0.033.  The  I)a  nf 
the  ui.  shea  red  liquid  decrease  hv  approximately  an  order  of  magnitude  on  increasing  the  density 
from  p  0.84-12  to  1.2ot>3.  However,  at  larger  shear  rates  the'  densitv  has  a  diminishing  influ- 
•  -nee  op  the  diffusion  coefficients.  The  high  shear  stresses  enhance  the  fluidity  and  create  a 
pud"  which  has  directionality.  All  the  diffusion  coefficient  components  rise  with  shear  rate 
.A  a  particular  density  because  the  structural  reorganization  that  takes  place  on  shearing  appar¬ 
ently  i  reutes  paths  along  which  particles  can  more  readily  move.  The  D  are  tvpicallv  two- 
thirds  <  f  the  Dx  and  Dz,  indicating  that  the  self-diffusion  is  favored  in  the  shearing  plane. 
'I  tins,  paradoxic  ally,  although  structural  aspects  of  the  the  licmids  under  shear  present  a  more 
solid-like  appearance,  other  more  dynamically  related  properties  such  as  self-diffusion  and 
shear  vise  osity  manifest  changes  which  are  associated  with  enhanced  fluidity. 

Thus  we  see  that  the  high  stresses  cause  extreme  distortions  of  the  system’s  structure. 
I  ne  mtermcdecular  forces  are  not  sufficiently  strong  to  support  the  enormous  stress  levels  that 
develop;  consequently  there  is  a  rupturing  and  reorganization  of  the  structure  resulting  jn 
enhanced  flow  behavior,  due  to  the  alignment  of  molecules  into  "glide'  planes  along  the  lino  of 
shear. 


These  structural  rearrangements  are  also  manifest  in  changes  in  the  so-called  normal 
pressure  components,  Pua  (-  -  Oaa)  Tvpicallv,  Pzz  >  Pxx  >  P  although  their  differenc  es  are 
small  when  compared  with  each  component's  change  on  shearing.  In  figure  A- 7  the  time  depen- 
dem  e  of  the  increase  in  P  ^  (t)  is  shown.  It  is  found  that  the  time  scale  of  the  development  el 
the  normal  pressure  components  is  longer  than  that  of  the'  shear  stress. 

file  difference  in  the  viscoelastic  and  structural  reorganization  times  can  he  better 
und.  (stood  b\  further  consideration  of  the  radial  distribution  functions. 

1  he  time  dependence  of  gac,  (r , t)  and  gag  !r,t)  reveal  that  structural  reorganization 
having  the  symmetry  of  shear,  <nz/t£'  >,  is  much  faster  than  that  of  symmetry.  Within 

;h  l.uf-r  serie-,  <x^/r^>  and  <z^/r^>  evolve  more  rapidlv  to  their  steady  state  values 
Hi  a  "yt  /r‘L  ’  For  example,  when  P  =  1.01304  and  e  =  0.7108,  the  former  three  averages  have 
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nearly  reached  their  steady  state  values  by  t  =  0.25  whereas  that  of  the  <y  /r  >  summation 
hardly  differs  from  zero  by  that  time.  Thus  we  find  that  when  a  liquid  is  sheared,  the  shearing 
forces  alter  its  structure  so  that  the  steady  state  stress  attained  corresponds  to  a  new  lower 
viscosity,  not  that  characterizing  the  starting  liquid.  At  high  rates  of  shear  the  structural 
evolution  necessary  to  attain  this  new  viscosity  is  significantly  slower  than  the  viscoelastic 
shear  relaxation.  Viewed  simply  this  means  that  the  original  viscosity  can  remain  for  a  suffi¬ 
ciently  long  time  to  enable  the  shear  stress  to  exceed  its  steady  state  value. 


Table  A  -2.  The  Density  and  Shear  Rate  Dependence  of  the  Directional  Diffusion  Coefficients, 

Da  ,  Obtained  from  the  MD  Experiments. 
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It  has  come  to  be  increasingly  recognized  that  the  process  of  fracture  involves  time 
dependent  nonlinear  viscoelastic  behavior.  In  order  to  better  understand  the  time  dependent 
nonlinear  viscoelasticity,  the  MD  calculations  were  conducted  with  a  time  varying  strain  rate. 
In  what  follows  we  describe  the  viscoelastic  MD  "experiments"  and  the  results.  What  is  most 
encouraging  is  that  in  every  case  the  behavior  of  real  liquids  is  mimicked  by  the  MD  results. 

We  consider  the  dynamical  response  of  the  LJ  material  to  an  imposed  shear  that  is 
"switched  on"  suddenly  using  the  Homogeneous  Shear  method.  Two  types  of  MD  experiments 
are  of  interest:  (a)  the  response  of  the  system  to  a  steady  shearing  rate  (uniform  velocity 
gradient)  «  =  d  v  /  d  z  that  is  initiated  at  time  t  =  0;  and  (b)  the  response  of  the  system  to  a 
constant  shear  strain  (a  pulse  of  shearing  rate)  imposed  at  time  t  =  t'.  Linear  response  theory 
provide,  us  with  relations  between  the  two  types  of  responses:  the  shear  stress  response  to  a 
step  of  shear  strain  can  be  written  as 


<7  (t)/«  =  (6) 

where  <l‘  (t)  is  the  normalized  (<f>(0)  =  1)  shear  stress  relaxation  function.  In  terms  of  these  pa¬ 
rameters  the  shear  stress  response  to  a  "small"  steady  shear  rate  imposed  at  t  =  0  is, 


a  (t)/e  =  G<dJ0  dt'  <t>  (  t  ’ 


(7) 


Since  the  steady  state  (t  ,  »)  value  of  this  is  just  the  viscosity  of  the  material,  tj  ,  we  have, 

no  -  G~/o  dt’  *  (t’>  =  G“T> 

The  last  equality  serves  to  define  the  shear  realization  time  r. 
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Experiments  of  the  type  described  as  (a)  above  were  then  carried  out  as  a  function  of 
t  for  values  of  £  well  beyond  the  range  for  which  linear  behavior  can  be  expected.  A  represen¬ 
tative  selection  of  the  results  is  displayed  in  figures  A-8  and  A-9  in  which  the  density  (con¬ 
stant  f.  1  and  shear  rate  (constant  density)  dependence  of  17  are  shown,  respectively. 

Several  features  of  these  figures  are  worthy  of  comment.  At  short  times  the  system 
behaves  largely  as  an  elastic  solid  in  that  the  stress  increases  linearly  with  time, 

o(t)  =  G  ^et.  (9) 


The  shear  rigidity  modulus  rises  linearly  with  pressure  at  roughly  the  same  rate  as  more  compli¬ 
cated  fluids,  as  is  revealed  in  figure  A-10.  Liquids  cannot  permanently  support  a  shear  strain 
and  consequently  the  stress  levels  off  to  a  value  determined  by  the  shear  viscosity  of  the 
sheared  state.  There  is  evident  "shear  thinning,"  or  a  decrease  of  viscosity  with  increasing 
dinar  rate.  The  predicted  value  of  the  stress  build-up  for  the  p  =  0.92862,  £  =  0.6963  state 
using  equation  (7),  shown  in  fig-are  A-8,  clearly  overestimates  the  actual  stress  obtained. 

Also,  a  maximum  in  the  stress  versus  time  curve  for  those  simulations  in  the  high 
c  nsitv/sh  ar  .-ate  regime  is  observed.  This  indicates  that  the  large  shear  deformations  are 
producing  structural  distortions  which  lag  in  time  somewhat  behind  the  growth  of  the  shear 
stress  hut  follow  rather  closely  the  observed  time  variation  of  the  normal  pressure  components 
m  the  materia!  This  relatively  slower  development  of  what  can  be  thought  of  as  transient 
stress-induced  glide  planes  is  consistent  with  the  observed  shear  thinning  and  stress  "overshoot" 
behavior  presented.  It  is  also  consistent  with  the  observation  that  the  stress  relaxation  function 
decays  more  rapidly  in  a  highly  sheared  system  than  in  an  equilibrium  system  (see  discussion  of 
g(  ^(t)  above).  A  comparison  of  the  two  types  of  behavior  is  presented  in  figure  A-ll.  The  form 
of  these  4,  i  t)  c  onsists  of  a  rapidly  (t  <  0.2)  decaying  portion  which  is  presumably  inertial  in 
origin  and  due  to  small  motions  of  the  molecules  within  the  random  network  structure.  Further 
stress  relaxation  requires  the  cooperative  rearrangement  of  larger  molecular  groups  and  is 
c  msequentlv  slower.  The  attenuation  of  d>(t)  with  increasing  shear  rate  is  indicative  of  the 
enhanced  stress  relieving  properties  of  the  structurally  reorganized  system.  The  stress  over¬ 
shoot  observed  in  nonlinear  viscoelastic  experiments  on  polymers  has  a  form  that  is  quite 
similar  to  the  curves  obtained  for  the  LJ  system.  While  the  microscopic  mechanisms  for  the 
two  systems  are  obviously  quite  different,  the  general  principle  that  associates  the  stress  over¬ 
shoot  phenomena  (as  well  as  shear  thinning)  with  a  time  dependent  structural  reorganization 
process  remains  quite  valid. 

In  this  context  it  is  perhaps  relevant  to  note  that  a  more  rapid  decay  of  the  stress 
relaxation  function  can  be  thought  of  as  both  a  suppression  of  the  long  relaxation  times 
associated  with  stress  relaxation  and  also  a  reduction  in  the  shear  rigidity  modulus,  at  low 
frequencies,  G,  (a  rather  well  known  phenomenon  in  polymeric  systems).* 

Thus  we  see  that  the  viscoelastic  studies  at  high  shear  rates  reflect  the  structural  reor¬ 
ganization  that  is  ultimately  responsible  for  failure  or  fracture  of  the  liquid.  The  fact  that 
every  property  "measured"  at  high  shear  rates  are  found  in  real  and  far  more  complicated 
liquids  suggests  strongly  that  MD  studies  cam  be  utilized  to  investigate  the  comminution  process 
in  those  liquids. 


*  V.  (!.  Yanovski  in  Relaxation  Phenomena  in  Polymers,  ed.  G.  Bartenev  and  Y.  V.  Zelenev 
(Khinhva,  Leningrad,  1972),  p.  113. 


Appendix 


33 


Experiments  of  the  type  described  as  (a)  above  were  then  carried  out  as  a  function  of 
t  for  values  of  £  well  beyond  the  range  for  which  linear  behavior  can  be  expected.  A  represen¬ 
tative  selection  of  the  results  is  displayed  in  figures  A-8  and  A-9  in  which  the  density  (con¬ 
stant  f.  1  and  shear  rate  (constant  density)  dependence  of  17  are  shown,  respectively. 

Several  features  of  these  figures  are  worthy  of  comment.  At  short  times  the  system 
behaves  largely  as  an  elastic  solid  in  that  the  stress  increases  linearly  with  time, 

o(t)  =  G  ^et.  (9) 


The  shear  rigidity  modulus  rises  linearly  with  pressure  at  roughly  the  same  rate  as  more  compli¬ 
cated  fluids,  as  is  revealed  in  figure  A-10.  Liquids  cannot  permanently  support  a  shear  strain 
and  consequently  the  stress  levels  off  to  a  value  determined  by  the  shear  viscosity  of  the 
sheared  state.  There  is  evident  "shear  thinning,"  or  a  decrease  of  viscosity  with  increasing 
dinar  rate.  The  predicted  value  of  the  stress  build-up  for  the  p  =  0.92862,  £  =  0.6963  state 
using  equation  (7),  shown  in  fig-are  A-8,  clearly  overestimates  the  actual  stress  obtained. 

Also,  a  maximum  in  the  stress  versus  time  curve  for  those  simulations  in  the  high 
c  nsitv/sh  ar  .-ate  regime  is  observed.  This  indicates  that  the  large  shear  deformations  are 
producing  structural  distortions  which  lag  in  time  somewhat  behind  the  growth  of  the  shear 
stress  hut  follow  rather  closely  the  observed  time  variation  of  the  normal  pressure  components 
m  the  materia!  This  relatively  slower  development  of  what  can  be  thought  of  as  transient 
stress-induced  glide  planes  is  consistent  with  the  observed  shear  thinning  and  stress  "overshoot" 
behavior  presented.  It  is  also  consistent  with  the  observation  that  the  stress  relaxation  function 
decays  more  rapidly  in  a  highly  sheared  system  than  in  an  equilibrium  system  (see  discussion  of 
g(  ^(t)  above).  A  comparison  of  the  two  types  of  behavior  is  presented  in  figure  A-ll.  The  form 
of  these  4,  i  t)  c  onsists  of  a  rapidly  (t  <  0.2)  decaying  portion  which  is  presumably  inertial  in 
origin  and  due  to  small  motions  of  the  molecules  within  the  random  network  structure.  Further 
stress  relaxation  requires  the  cooperative  rearrangement  of  larger  molecular  groups  and  is 
c  msequentlv  slower.  The  attenuation  of  d>(t)  with  increasing  shear  rate  is  indicative  of  the 
enhanced  stress  relieving  properties  of  the  structurally  reorganized  system.  The  stress  over¬ 
shoot  observed  in  nonlinear  viscoelastic  experiments  on  polymers  has  a  form  that  is  quite 
similar  to  the  curves  obtained  for  the  LJ  system.  While  the  microscopic  mechanisms  for  the 
two  systems  are  obviously  quite  different,  the  general  principle  that  associates  the  stress  over¬ 
shoot  phenomena  (as  well  as  shear  thinning)  with  a  time  dependent  structural  reorganization 
process  remains  quite  valid. 

In  this  context  it  is  perhaps  relevant  to  note  that  a  more  rapid  decay  of  the  stress 
relaxation  function  can  be  thought  of  as  both  a  suppression  of  the  long  relaxation  times 
associated  with  stress  relaxation  and  also  a  reduction  in  the  shear  rigidity  modulus,  at  low 
frequencies,  G,  (a  rather  well  known  phenomenon  in  polymeric  systems).* 

Thus  we  see  that  the  viscoelastic  studies  at  high  shear  rates  reflect  the  structural  reor¬ 
ganization  that  is  ultimately  responsible  for  failure  or  fracture  of  the  liquid.  The  fact  that 
every  property  "measured"  at  high  shear  rates  are  found  in  real  and  far  more  complicated 
liquids  suggests  strongly  that  MD  studies  cam  be  utilized  to  investigate  the  comminution  process 
in  those  liquids. 


*  V.  (!.  Yanovski  in  Relaxation  Phenomena  in  Polymers,  ed.  G.  Bartenev  and  Y.  V.  Zelenev 
(Khinhva,  Leningrad,  1972),  p.  113. 


Appendix 


33 


no*i 


aftl 

o 


biyure  4-8. 


Appendix 


t 

The  Time  Dependence  of  the  Reduced  Stress,  a(t)/e: _ p  =  0.8442,  eo  =  I 

.  (a)  Actual  .  (b)  Predicted  Using  Equation  (4)  p  =  0.92862, 

0.6963; _ p  =  1.01304,  kQ  =  0.7168;  __P  =  1.2663,  eo  =  0.7721. 
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Figure  A-9.  The  Time  Dependence  of  the  Reduced  Stress,  CJ(t)/Gg,  for  the  p 
State: _ eQ  =  0.0896, _ eQ  =0.3584,  ....  e  =  0.7168. 


=  1.01304 


Append  i  x 


35 


po.|q 


20 


t 


Figure  A-9.  The  Time  Dependence  of  the  Reduced  Stress,  CJ(t)/Gg,  for  the  p 
State: _ eQ  =  0.0896, _ eQ  =0.3584,  ....  e  =  0.7168. 


=  1.01304 


Append  i  x 


35 


The  Pressure  Dependence  of  G#:  ASantotrac  40,  O  Di(2-ethylhexyl)  Phthalate 
from  a  Twin-Disk  Apparatus.  Solid  line  =  (0.293  +  1.70P)  GPa  di(2- 

e)phthalate.  The  filled  in  O  is  obtained  from  MHS  MD. 
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from  a  Twin-Disk  Apparatus.  Solid  line  =  (0.293  +  1.70P)  GPa  di(2- 

e)phthalate.  The  filled  in  O  is  obtained  from  MHS  MD. 
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Commandant 
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Commandant 

US  Army  Missile  A  Munitions  Center 
and  School 
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Commandant 

USAMP&CS/TCAFM 
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Commander 
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Director,  Development  Center 
Marine  Corps  Development  and 
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ATTN:  Fire  Power  Division  1 

Quantico,  VA  22134 


41 


US  ARMY  ARMAMENT  MATERIEL  READINESS 
COMMAND 


Commander 

US  Army  Armament  Materiel 

ATTN:  DRSAR-ASN 

ATTN:  DRSAR-SF 

ATTN:  DRSAR-SR 

Rock.  Island,  IL  61299 


Readiness  Command 
1 
1 
1 


Commander 
USA  ARRCOM 

ATTN:  SARTE  1 

Aberdeen  Proving  Ground,  MD  21010 


Comma  n  der 

US  Army  Dugway  Proving  Ground 

ATTN:  Technical  Library  (Docu  Sect)  1 

Ougway,  UT  8'*922 


US  ARMY  TRAINING  A  DOCTRINE  COMMAND 


Commandant 

US  Army  Infantry  School 

ATTN:  NBC  Division  1 

Fort  Benning,  GA  31905 

Commandant 

US  Army  Missile  A  Munitions  Center 
and  School 

ATTN:  ATSK-OT-MU-EOD  1 

Restone  Arsenal,  AL  35809 

Commandant 

USAMP&CS/TCAFM 

ATTN:  AT  2N-CM-CDM  1 

Fort  McClellan,  AL  36205 

Commander 

US  Army  Infantry  Center 

ATTN:  ATSH-CD-MS-C  1 

Fort  Benning,  GA  31905 

Commander 

US  Army  Infantry  Center 
Directorate  of  Plans  A  Training 
ATTN:  ATZB-DPT-PO-NBC  1 

Fort  Benning,  GA  31905 

Comma  nder 

USA  Training  and  Doctrine  Command 

ATTN:  ATCD-Z  I 

Fort  Monroe,  VA  23651 


Commander 

USA  Combined  Arms  Center  and 
Fort  Leavenworth 

ATTN:  ATZL-CA-COG  1 

ATTN:  ATZL-CAM-IM  I 

Fort  Leavenworth,  KS  66027 

Commander 

US  Army  TRADOC  System  Analysis  Activity 
ATTN:  ATAA-SL  1 

White  Sands  Missile  Range,  NM  88002 

US  ARMY  TEST  A  EVALUATION  COMMAND 

Commander 

US  Army  Test  A  Evaluation  Command 

ATTN:  DRSTE-CT-T  1 

Aberdeen  Proving  Ground,  MD  21005 

DEPARTMENT  OF  THE  NAVY 

Chief  of  Naval  Research 

ATTN:  Code  443  1 

800  N.  Quincy  Street 
Arlington,  VA  22217 

Commander 

Naval  Explosive  Ordnance  Disposal  Facility 
ATTN:  Army  Chemical  Officer  (Code  AC-3)  1 

Indian  Head,  MD  20640 

Commander 

Naval  Surface  Weapons  Center 

Code  G51  1 

Dahlgren,  VA  22448 

Chief,  Bureau  of  Medicine  A  Surgery 
Department  of  the  Navy 

ATTN:  MED  3C33  1 

Washington,  DC  20372 

Commander 

Nava  I  Weapons  Center 

ATTN:  Technical  Library  (Code  343)  1 

China  Lake,  CA  93555 

US  MAR  I NE  CORPS 

Director,  Development  Center 
Marine  Corps  Development  and 
Education  Command 

ATTN:  Fire  Power  Division  1 

Quantico,  VA  22134 


41 


department  of  the  air  force 

hQ  Foreign  Technology  Division  (AFSC) 

ATTN:  TQTR  1 

Wr  i  glit-Patterson  AFB,  OH  45433 


Commander 

Aeronautical  Systems  Division 
ATlN:  ASD/AELC  1 

ATTN:  ASD/AESD  1 

Wr I ght-Patterson  AFB,  OH  45433 

HQ  AFIC/LOWMM  1 

Wr i ght- Patterson  AFB,  OH  45433 

HQ,  aFSC/'DNL  1 

Andrews  AFB,  MD  20334 

HQ  AMO/RD 

ATTN:  Chemical  Defense  OPR  1 

Brooks  AFB,  TX  78235 

N0RA0  Combat  Operations  Center 

ATTN:  DOUN  I 

Cheyenne  Mtn  Complex,  CO  80914 


Air  Force  Aerospace  Medical  Research 
La  boratory 

ATTN:  AFAMRL/HE  (  Or.  C.R.  Replogle)  1 

Wr i ght-Patterson  AFB,  OH  45433 

HQ  AFTEC/SGB  1 

Kirtland  AFB,  NM  87117 

OUTSIDE  AGENCIES 

Battel le,  Columbus  Laboratories 
ATTN:  TACTEC  I 

505  Kl ng  Avenue 
Columbus,  OH  >201 

Toxicology  Information  Center,  WG  1008 
National  Research  Council  1 

2101  Constitution  Ave. ,  NW 
Washington,  DC  20418 

US  Public  Health  Service 
Center  for  Disease  Control 

A1TN;  Lewis  Webb,  Jr.  1 

Building  4,  Room  232 
Atlanta,  GA  30333 
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